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We investigate cross sections for neutrino- 12 C exclusive scattering and for muon capture on 12 C 
using wave functions obtained in the ab initio no-core shell model. In our parameter-free calcu- 
lations with basis spaces up to the 6hQ we show that realistic nucleon-nucleon interactions, like 
e.g. the CD-Bonn, under predict the experimental cross sections by more than a factor of two. 
By including a realistic three-body interaction, Tucson-Melbourne TM'(99), the cross sections are 
enhanced significantly and a much better agreement with experiment is achieved. At the same time, 
the TM'(99) interaction improves the calculated level ordering in 12 C. The comparison between the 
CD-Bonn and the three-body calculations provides strong confirmation for the need to include a 
realistic three-body interaction to account for the spin-orbit strength in p-shell nuclei. 

PACS numbers: 25.30.Pt, 21.60.Cs, 21.30.Fe, 27.20. +n 

The Gamow- Teller (GT) transition from the ground state of 12 C to the 1+ T=l isobar triplet ( 12 B g . s ., 12 C(15.11 
MeV), 12 N g . s .) is a very sensitive test of nuclear structure models for mass 12 and, particularly, of the strength of the 
spin-orbit interaction. The two most common p-shell approximations for the structure of the ground state of 12 C ((a) 
the p-shell equivalent of a L=0 S=0 three alpha-cluster structure and (b) the closed p 3 / 2 shell structure) give very 
different (indeed opposite) predictions for the B(GT) strength to T=l 1 + triplet. In the p-shell alpha-cluster limit the 
ground state of carbon has good SU(4) symmetry [444] and the Gamow- Teller transition is forbidden because there 
does not exist a 1 + T=l state with [444] symmetry and <tt operator cannot change SU(4) symmetry. This translates 
into an exact cancellation between the different P1/2 and p 3 / 2 transition amplitudes. The observed transition strength 
requires the inclusion of higher SU(4) components in the wave functions and the breaking of the cancellation is quite 
sensitive to the assumed spin-orbit interaction. In the the jj-coupling limit, where one assumes that the ground state 
of 12 C is described by a closed p 3 / 2 shell, the transition to the T=l 1 + state is pure p 3 / 2 — * Pi/2- No cancellations 
between different transition amplitudes are allowed and the transition strength is over estimated by almost a factor 
of 6. When RPA correlations are included in the initial and final states the situation improves somewhat, but the 
transition remains over-estimated by about a factor of 4 0, |2| . The strong contrast between the predictions of the 
pure jj-coupling and the pure SU(4) limits makes this Gamow- Teller transition an ideal test case for the strength of 
the spin-orbit interaction and for model wave functions of mass 12. 

In this letter we present the predictions of no-core shell model (NCSM) |3| calculations of 12 C for the T=l 1 + 
transition in 12 C. We examine inelastic electron scattering to the 15.11 MeV state of 12 C, muon capture to the 
ground state of 12 B, and neutrino scattering to the ground state of 12 N. These different electroweak reactions probe 
different momentum transfers and comparisons between theory and experiment allow us to test the convergence of 
the no-core shell model with increasingbasis size, up to 6h£l. We also investigate the contributions of a three-nucleon 
force since it is now well established 0, 0, 0] that realistic nucleon-nucleon interactions alone account for only half 
the observed p-shell splitting, while the rest arises form two-pion exchange between three or more nucleons. In the 
present calculations we include a realistic chiral-symmetry-based three-nucleon interaction (TNI) , Tucson-Melbourne 
TM'(99)0. 

A detailed description of the NCSM approach was presented, e.g. in Refs. 0,13. Here, we simply present extensions 
and modifications needed when a genuine TNI is included. The starting Hamiltonian is Ha — -j ^2i<j ^'2!^ — 
^2t<j ^NN.ij + ^2f<j<k ^NNN,r/fc, where Vnn.ij is the nucleon-nucleon (NN) interaction and VNNN.ijfc is the TNI. We 
employ a large but finite harmonic-oscillator (HO) basis. Due to properties of the realistic nuclear interaction we 
have to derive an effective interaction appropriate for the selected finite basis space. To facilitate this, we modify 
the Hamiltonian by adding to it the center-of-mass (CM) HO Hamiltonian Fcm = 7cm + Uqm, where Ucm — 
iAmfl 2 R 2 , R = 4 Ylt=i ^i- The effect of the HO CM Hamiltonian will later be subtracted in the final many-body 
calculation. The modified Hamiltonian can be cast into the form = Ha + -Hcm = Sj=i hi + ^>t<j ^ij ' A + 

J2f<j<k ^NNN,»ifc, where h t = ^ + ±mf2 2 rf and V^' A = Vnn,^ - 2 §x-(^ - ^) 2 - Next w e divide the A-nucleon 
infinite HO basis space into the finite active space (P) comprising all states up to A^ max HO excitations above the 
unperturbed ground state and the excluded spaces (Q = 1 — P). The basic idea of the NCSM approach is to apply a 
unitary transformation on the modified Hamiltonian, e^ s H^e s such that Qe~ s H^e s P = 0. If such a transformation 
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is found, the effective Hamiltonian that exactly reproduces a subset of eigenstates of the full space Hamiltonian is 
given by H e g = Pe~ s H^e 3 P. This effective Hamiltonian contains up to A-body terms and it is essentially as difficult 
to construct it as to solve the full problem. Therefore, we apply this approach with a cluster approximation. When 
a genuine TNI is considered, the simplest cluster approximation produces a three-body effective interaction. The 
NCSM calculation is then performed in four steps: 

(i) We solve a three- nucleon system for all possible three-nucleon channels with the Hamiltonian H^,i.e., using 
hi + h 2 + h 3 + V 12 ' + V 13 ' + V 2 ^' + V^nnn,123- Consequently, the three nucleons feel a pseudo-mean field of the 
spectator nucleons generated by the HO CM potential. It is necessary to separate the three-body effective interaction 
contributions from the TNI and from the two-nucleon interaction. Therefore, we need to find three-nucleon solutions 
for the Hamiltonian with and without the Vnnn 123 TNI term. The three-nucleon solutions are obtained by procedures 
described in Refs. Q (without TNI) and [l(| (with TNI). 

(ii) We construct the unitary transformation corres pon ding to the choice of the active basis space P from the 
three-nucleon solutions using the Lee-Suzuki procedure [llL Ha ] . The three-body effective interaction is then obtained 
as ^3cff,l 2 N 3 NN = P[e- s ™+ N ™(/ il + h 2 + h 3 + V$ A + V^ A + V 2 Y + y N NN,i23)e s ™+™ N - fa + h 2 + h 3 )]P and 
Kjeff.123 = P[e- s "™(/ii + h 2 + h 3 + VH' A + V^' A + V 2 n 3 ' A )e Sm - {hi + h 2 + h 3 )]P. The three-body effective interaction 
contribution from the TNI is then defined as V 3 N JJ 23 = V^j^ NN - V 3 N J 123 . 

(iii) As the three-body effective interactions are derived in the Jacobi-coordinate HO basis but the A = 12 cal- 
culations will be performed in a Cartesian-coordinate single-particle Slater-determinant m-scheme basis, we need 
to perform a suitable transformation of the interactions. This transformation is a generalization of the well-known 
transformation on the two-body level that depends on HO Brody-Moshinsky brackets. 

(iv) We solve the Schrodinger equation for the A = 12 nucleon system using the Hamiltonian H^ cS — Y*,t=i + 

A^2 J2t<j<k Vitff.ijk + J2t<j<k ^3eff,yfc' where the factor takes care of overcounting the contribution from the 
two-nucleon interaction. At this point we also subtract the i?CM- The A = 12 nucleon calculation is then performed 
using the Many-Fermion Dynamics shell model code ^| generalized to handle three-body interactions. Eventually, 
the transition densities are computed that serve as an input for evaluating our selected observables. 

Detailed 12 C NCSM calculations using realistic two-nucleon interactions were reported in Ref. Here we extend 
those calculations by including the TNI and reach the Ahtt (6h£l) basis in calculations with (without) the TNI. 
In Table [I] we summarize some of our results. In general, in addition to increase of binding energy, we observe a 
substantial sensitivity of the low-lying spectra to the presence of the TNI and a trend toward level-ordering and 
level-spacing improvement in comparison to experiment. The sensitivity is the largest for states where the spin-orbit 
interaction strength is known to play a role. Note the correct ordering of the 1+0 <-> 4+0 states and ordering and 
spacing improvement of the lowest T=l states. 

The significant increase in the spin-orbit splitting obtained from the inclusion of the TNI is seen most strikingly 
in the predicted cross sections to the T=l 1 + states. In all our electron scattering and weak interactions results 
here, only one-body currents are included and the bare operators are used. Dubach and Haxton |17| have shown that 
at high-momentum-transfers it is necessary to include two-body meson- exchange currents to describe the transverse 
magnetic electron scattering form factor for excitation of the 15.11 MeV state. However, a reasonable description of 
the form factor up to momentum transfers of about 200 MeV/c can be obtained with a one-body current. 

Figurenshows a comparison between the form factors predicted by the NCSM and experiment. The experimental 
data are represented by the black circles, which represent a fit to the data assuming only a one-body current obtained 
by Dubach and Haxton [l^. The theoretical curves shown are NCSM results for 2M1, AMI, 6Ml using the CD-Bonn 
NN interaction and a AMI calculation using the AV8' plus the TM'(99) realistic TNI. The qualitative features of our 
results are seen by looking at the height of the first maximum and the position of the minimum. With two-body 
interactions alone, the change in transition form factor from AhQ to 6h£l is small compared to the differences between 
theory and experiment. The magnitude of the form factor is too low and minimum occurs too far out in momentum. 
When the TNI is included a significant improvement is seen in both the shape and magnitude of the theoretical form 
factor. The magnitude of the form factor up to the first maximum is close but somewhat lower than experiment. The 
shape of the form factor is also improved but it is still stretched out too far in momentum. Comparing the AMI 3-body 
calculation with the AMI CD-Bonn calculation the magnitude of the form factor at the peak has increased by about 
75%, and the position of the first minimum has shifted from q m i n ~ 400 MeV/c to 360 MeV/c. This improvement is 
almost entirely due to the improved strength of the spin-orbit splitting when the TNI is included. 

The 15.11 MeV state was included in the fit to the p-shell interaction by Cohen and Kurath [l8j (CK), and the CK 
interaction probably represents the best description of this transition using a globally fitted p-shell interaction. Our 
AV8'+TM'(99) form factor agrees well with the CK prediction, when for consistency between the two calculations, 
we use b=1.663 fm. 

The conclusion drawn from the transverse magnetic form factor results is further supported by our B(M1; 1 + 1 — > 
+ 0) results presented in Table [I]and Fig. [3 The calculations with 2-body forces show saturation and underpredict 
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the experiment by almost a factor of three. By including the TNI, the B(M1) value increases dramatically. We fully 
expect that further increases in the basis size will produce results with TNI close to experiment. For smaller basis 
sizes, effective transition operators may be important and work in this direction is underway. 

Table |n] shows the comparison between the theoretical and experimental neutrino scattering cross section for the 
same selections of Hamiltonians and basis spaces. These results show a similar trend to the electron scattering results 
above. In this case, the neutrino spectrum for electron neutrinos from decay-at-rest (DAR) of the pion peaks around 
30 MeV and the average momentum transfer is about 40 MeV/c. The CD-Bonn interaction (without TNI) results 
indicate an approach to convergence by 6hQ but experiment is under-predicted by about a factor of 2.4. When the 
TNI is included with the AV8' interaction the predicted cross section is only 30% lower than experiment. Based on 
the similarity of trends with the electron scattering results, we anticipate that when the model space is eventually 
expanded to 6hil theory would be within 15% of experiment. The substitution of AV8' for CD-Bonn in the calculations 
with TNI is expected to be of minor consequence. 

Muon capture involves a higher momentum transfer than the (v e , e~) reaction and the average momentum transfer 
is q ~ 100 MeV/c. By QUO, the CD-Bonn calculations show signs of converging yet experiment is underestimated 
a factor of 2.6. The inclusion of the TNI shows a significant improvement and, for 4fi,f2, theory is 34% lower than 
experiment. Again extrapolating using the trends of the inelastic electron scattering results suggests that a QMl 
calculation that included a realistic TNI would come within 20% of experiment. 

The (l/p, /i~) neutrino cross section to 12 N ff s . corresponds to the LSND muon neutrinos from decay-in-flight (DIF) 
of the pion. This spectrum involves neutrinos up to about 250 MeV, with a average neutrino energy of about 150 
MeV and an average momentum transfer of about 200 MeV/c. In this case the QhQ. CD-Bonn calculation is off by a 
factor of 1.8 compared with experiment. The AMI calculation that includes the 3-body TM'(99) interaction is, in fact, 
in agreement with experiment. However, based on the trends established above, this suggests that a larger model 
space may over-predict experiment. Examining the elastic scattering form factor suggests that the problem lies in the 
fact that at 200 MeV/c the predicted form factor is too large. Of course, as the model space is increased we expect 
the form factor to be shifted down in momentum. 

In conclusion, the transition from the 12 C S . S . to the T=l 1 + states in mass 12 is very sensitive to the strength 
of the spin orbit interaction. We have investigated neutrino- 12 C exclusive cross sections and muon capture on C 
as well as inelastic electron scattering using wave functions obtained in the ab initio NCSM. In our parameter-free 
calculations with basis spaces up to 6h£l we show that realistic NN interactions under predict the experimental weak 
interaction cross sections by more than a factor of two. At high momentum transfers around q ~ 200 MeV/c the 
electron scattering form factor is over-predicted and the position of the predicted minimum to close to 400 MeV/c, 
compared to the experimental minimum at q ~ 260 MeV/c. By including a realistic TNI the weak interaction cross 
sections are enhanced significantly which considerably improves agreement with experiment. The shape of the electron 
scattering form factor is also significantly improved, but the predicted form factor still peaks at too large a momentum 
transfer and is too large at the momentum transfers relevant to the LSND DIF cross section. The difference between 
the observed and predicted shape for the (e,e') form factors and the very different momentum transfers involved in the 
three weak processes examined here imply that a single experiment /theory scale factor cannot be defined for all three. 
The comparison between the CD-Bonn and the three-body calculations discussed here provide a strong confirmation 
of the need to include a realistic three-body interaction to account for the spin-orbit strength in p-shell nuclei. 

This work was performed in part under the auspices of the U. S. Department of Energy by the University of Cali- 
fornia, Lawrence Livermore National Laboratory under contract No. W-7405-Eng-48. This work was also supported 
in part by USDOE DE-FG02-87ER40371. 
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FIG. 2: B(M1) values, in fi%, of the 0+0 -> 1+1 transition in 12 C. For details see Table[j] 
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AV8'+TM'(99) 
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CD-Bonn 






basis space 




4htl 


4hti 


6hQ, 


4hSl 


2UQ. 


OhQ. 


\E gs \ [MeV] 


92.162 


91.963 


85.944 


85.630 


88.518 


92.375 


104.947 


Q 2+ [e fin 2 ] 


+6(3) 


4.288 


4.613 


4.717 


4.532 


4.430 


4.253 


£ x (2+0) [MeV] 


4.439 


3.603 


3.427 


3.612 


3.697 


3.837 


3.734 


£ x (l+0) [MeV] 


12.710 


11.280 


13.926 


13.930 


14.140 


14.524 


13.866 


£ x (4+0) [MeV] 


14.083 


13.517 


12.272 


13.110 


13.356 


13.638 


12.406 


S x (l+1) [MeV] 


15.110 


16.221 


16.364 


16.064 


16.165 


16.291 


15.290 


E K (2+1) [MeV] 


16.106 


16.467 


17.712 


17.409 


17.717 


17.945 


15.970 


£ x (0+l) [MeV] 


17.760 


17.116 


16.213 


16.534 


16.619 


16.493 


14.698 


B(E2;2+0 -> 0+0) 


7.59(42) 


4.146 


4.765 


5.019 


4.624 


4.412 


4.092 


B(M1;1+1 -> 0+0) 


0.951(20) 


0.645 


0.305 


0.384 


0.355 


0.280 


0.158 


B(E2;2+1 -> 0+0) 


0.65(13) 


0.430 


0.247 


0.309 


0.283 


0.015 


0.002 



TABLE I: Experimental and calculated properties of 12 C. The units are e 2 fm 4 (fi%) for B(E2) (B(M1)). Three-body effective 
interactions derived from the AV8' [T3 ] and AV8'+TM'(99) and two-body effective interactions derived from the CD-Bonn [l5| 
NN potential and a HO frequency of fiQ. — 15 MeV were used. The TM'(99) parameters are given in Refs. [?1 lir| with the 
cutoff set to A = 4.7. The experimental values are from Ref. [T3| . By extrapolating our results we predict the CD-Bonn 12 C 
binding energy to be ~ 80 ± 2 MeV. 



Interaction 


CD-Bonn 2hQ, CD-Bonn 4hQ, CD-Bonn 6hQ, AV8'+TM'(99) ATiQ, 


experiment 


("mi M~) 
/i-capture 


2.27 3.2 3.69 6.8 
0.168 0.275 0.312 0.537 
1.46 2.07 2.38 4.43 


8.9±0.3±0.9 [19J 
0.56±0.08±0.1 [20] 
6.0±0.4 [21] 



TABLE II: Predicted weak interaction rates for the 12 C^ T=l 1 + transitions. The units are 10 42 cm 2 for the (v e ,e ) DAR 
cross section, 10 _40 cm 2 for (v/j,,fj,~) DIF cross section and 10 3 sec _1 for muon capture. 



